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Abstract 

We study a model due to J.M. Lasry and P.L. Lions describing the evolution of a 
scalar price which is realized as a free boundary in a ID diffusion equation with dynam- 
ically evolving, non-standard sources. We establish global existence and uniqueness. 



Introduction 



Here we are concerned with the following PDE: 

ft - fxx = [5p(t)+a - '5p(t)_a] fx{p{t), t) in (-1, 1) X [0, co); 

(P) < /x(l,t)=/x(-l,<) = 0. 



[ /(x,0) = //(z). 

where p{t) = {x : f{x, t) = 0} presumed, for a.e. t, to be a singleton, and a ~ 
min{a, \p{t) ± 1|} with a < 1. 

The model with a = a was invented in [5] and, as explained therein (see also [S],[7]) is 
purported to describe the dynamic evolution of a price p(t) as influenced by a population 
of buyers and sellers. In this initial reference, the existence of solutions was discussed, 
mostly in the context of a non-compact domain. 

While the model on compact domains was featured (strictly speaking, the model on 
R does not make economic sense) there was no proviso for the circumstance | l±p(t) | < a. 
Our modification using a provides this definition and later (see Lemma 12. 2p allows us 
to ensure that p(t) stays away from the domain boundaries at .t = ±1. In terms of the 
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model, our modification can be viewed as a "rescue plan" to prevent prices from severe 
deflation or inflation. 

Recently, [S] the problem was solved completely for the case of symmetric initial 
data and in the work [7], global existence, uniqueness and stability was established 
for initial data sufficiently close (in a certain sense) to the piccewise linear equilibrium 
solution. Finally, contemporaneous to the present work, a regularized version of (P) 
for the non-compact case is reinvestigated in [13| . A complete derivation of uniqueness 
for short times - roughly the equivalent of our §1 - is presented therein. 

Notwithstanding the benign appearance of (P), the system contains intrinsic and 
convoluted non-linearities. Indeed, the driving term at the sources is the gradient 
at the dynamically generated zero (free boundary) of / - which in turn controls the 
location of the sources. Thus, a central technical issue is to establish non-degeneracy 
at the free boundary and thereby some degree of control for its motion. E.g., in this 
context Hopf 's lemma, while useful, is not immediately decisive without some additional 
regularity information at the boundary. 

We consider the initial data // e C^([— 1, 1]) satisfying the following: 

(i) {fi{x) = 0} = {pi}, fi{x) > for X < PI and fi{x) < for x > p/. 

(ii) a,//(-i) = a,/K+i) = o. 

(iii) Given Xj := —dxfi{pi) > 0, we must have 

-dxfi{x)>^Xi in {pj - ao,pj + ao) (0.1) 

for some < oq ^ a/4. 

It is worthwhile to notice that problem (P) satisfies the following conservation iden- 
tities 

p(t) ppi nl nl 

fix, t)dx^ / fi{x) dx = Mb > 0, / fix, t) dx = / fjix) dx ^ -Mp < 0. 
-1 J-i Jp{t) J PI 

(0.2) 

Note also that the free boundary moves with a velocity given by the formula: 

fxx{p{t),t) 
fxipit),t) 

and the flux across the free boundary is given by 

A(<) ■.= -d,fipit),t). 

The main result in this paper is written below: 



^It is noted - but not proven - that existence might be established via a connection to a stochastic 
interacting particle model. The utility of this connection is under investigation by the authors, particularly 
with regards the question of global stability. 
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Theorem 0.1. [Global existence and uniqueness of classical solution] 

Consider the system described on (P) with initial data fj satisfying the conditions 
(i)~(iii) above, and let us define := —{ — 1,1) x{0,oo). Then there is a unique function 
f{x,t) G L°°(fl) satisfying the following: 

(A) For any < t <T there exists r = r[T) with r > such that 



(a) f IS C°° in {{x,t) : \x - p{t)\ < r}; 

(b) X{t) > r; 

(c) p{t) e (-l + r,l-r). 

(B) f solves the first two equation of (P) in the classical sense (in terms of Duhamel's 
formula ) in fi, and f{x,t) uniformly converges to /(x,0) as t 0. 



1 Short times 

The preliminary results are based on the short-time contraction principle of the follow- 
ing iteration: given /„(.T,t) and Pn{t) such that Pn{t) '■= {x : fn{-,t) = 0} consists of a 
unique point for each i > 0, consider the function A„(t) defined by 

A„W := -ifn)APn{t),t). 

Let f{x,t) solve 



df 9^/ 

'dt ^ 'd^ ^ '^"^^^ ['5x=p„(t)-a - '^a;=p„(t)+a] > (I'l) 

= /.(l,t) = 0, (1.2) 
/(x,0) = fi{x). 

where in the above, it is tacitly assumed that for all t in (O,to)i Pn{t) is a single point. 
Then the solution to the above becomes fn+i{x,t), i.e., serves to define Pn+iit) and 

Specifically, we define the map 

$: L°°((0,io);X) ^ i°=((0,to);^), (1.3) 

fn fn+1, 

where X :— L°°{pj — a^^pj + oq) for some oq sufficiently small. 

First let us write the solution f{x,t) of p.2p using the Duhamel formula: 
.+1 

f{x,t) = / Tix,x';t)fj{x')dx' 

/■* (1 4) 

+ [T{x,{pn{t')-a);t-t')-r{x,{pn{t') + a);t-t')]Xn{t') dt' 
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where F denotes the fundamental solution appropriate for the domain (— 1,+1) with 
Neumann boundary conditions: 

oo 

r{x,x';t)^ ^ K{x^{2k+[-lpx'),t) (1.5) 

k— — oo 

with 



K{x,t) = 




The main result of this section is stated in the following: 
Theorem 1.1. There exists a time to depending only on 0'''^d Xj such that 

sup ||/-5llx<^ sup |l/„-5„|lx, 
te[o.to] ^ telo.to] 

where X :~ L°° {pi ~ ao,pi + oq) and the functions /, g solve il.S^) with right-hand side, 
respectively, {f„,p„) and (5„,g„). 



Corollary 1.2. There exists to > depending only on jl//jlLoc(_x.i) "^"c? A/ such that 
(P) has a solution f,p for all times t G [0, ig]- Moreover, f is smooth accross the free 
boundary, and 

Xito) > ^ > 0. 



The proof of theorem II. II is a consequence of the following series of results: 
Lemma 1.3. For all x and y such that \x ^ y\ < a/4, it holds that 

T{x,y±a,t)+ \Ta;{x,y±a,t)\+ \r^x{x,y ± a,t)\ < G{a), 
where the constant G{a) depends only on a. 

Proof. It follows from elementary analytical considerations. We show it for T^, the rest 
are similar. Writing 

x-{2k+[-\\\^\x') 

tfh = , 

2Vt^ 

we have 

with Gi a constant. Obviously, for any fc, we can bound Gi0^c~^'= < G2 for another 
constant G2, whatever the value of 'dk might be. 

Next, since for the relevant x' = y±a we have that > a/2 and — 2| > g^a 

for a constant ^3, we may sum the series replacing a; ± x' by the relevant worst case 
scenarios. It is concluded that jr^;! - unintegrated in the t-variablc - is bounded by a 
finite constant (which depends on a, but not on x, x' or t). □ 
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Lemma 1.4. Consider 



fx{x,t) = I T,{x,x'-t)fi[x')dx' 



+ / \rx{x,pn{t') - a]t-t') -Tx{x,pn{t') ~\- a;t-t')\\n{t') dt' 
Jo 



If t is a time such that \pi — Pn{t)\ < cLq < a/4 for all t <t, then we can estimate 

dhix, t) 



dx 



<G{a) / \\n{s)\ds, 





for all t < t and x € (p/ — clq,pi + gq), where A/ > 0, pi are as described in the 
hypothesis (i) - (Hi), and G{a) is a constant depending only on a. 

Moreover, for x G {pi ~ ao,pi + ao), the linear term can be estimated as 

dli(x, t) 1 , t , , 

\ < -T^Aje-*, forallt<^. 
dx 2 



Proof. The second estimates is a direct consequence of Harnack inequality. Indeed, 
— ^ is a solution of an homogeneous heat equation with initial data —dxfiix), that 
has a definitive sign (and lower bound) throughout {pi — ao,p/ + ao). Thence we can 
apply the classical parabolic Harnack inequality from [T3] to — 



The first estimate follows easily from Lemma ll.3[ just taking into account that 
\x-Pn{t)\ < \x -pi\ + \pi -Pn{t)\ < 2ao << a. 



□ 



Corollary 1.5. Suppose, in the n stage of the iteration, that for all t < t we have 

|A„(0| <2M^zH) and b/ - p„(i)| < ao (1.7) 

Then we can find the following bound for the next step f = fn+i, given in p.4p ; there 
exists a time to < t, depending on \\fi\\i^ao(^_i j^-j and Xj , such that 

M[~f,ix,t)]>^, 
t,x 4 

for all t <tQ and \x — pi\ < miiijao, ao}. 
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Proof. This is an obvious consequence of the two bounds from Lemma 11.41 and the 
hypothesis ()1.7p . Indeed we may write, uniformly for x G {pi — ao.pi + ao) 

(/„+i),(x,t) <-iA/e-* + 4G(a)||MUoo(_i,i)Vt. 

It is clear that for t < with tg sufficiently small, and depending only on the quotient 
II//|Ilcxj(_i 1) /A/, the desired bound will hold. It is noted that the time to docs not 
depend on the iteration coefficient n. □ 

In the next two lemmas we show a L°°-bound for fx and fxx in a neighborhood of 
the free boundary which is the required input for Corollary 1 1.51 

Lemma 1.6. Suppose that i is a time such that \pi — Pn{t)\ < for all t < t and that 
sup|(/„)^(x,t)| <2 "-^'"^°°(-^'^\ 0<t<t, \x-pi\<ao, (1.8) 



the 



Vt 



sup|(/„+i),(a:,t)| <2M^i^9^^, 0<t<t, \x-pi\< do. 



Proof. As before, we let 



fx{x,t)= / Tx{x,x']t)fi{x') dx 



+ / [Tx{x,{pn{t')~a)-t-t')-Tx{x,Pn{t') + a-t~t')]\n{t')dt' 
Jo 

dli{x,t) ^ cfis . 
dx dx ' 

Regularity estimates for caloric functions (see [1], Chapter V, Theorem 8.1) imply 
dli{x, t) 



sup 



dx 



< -jz , for < oq. 



For the nonlinear part we make use of the estimates for the Green's function given in 
Lemma [T751 and the hypothesis p.9p . It holds 



sup 



dhix, t) 



dx 



<G{a) |A„(s)| ds<4G(a)||MUoo(_i,i)Vt, \x-pi\<do. 



As soon as t is small enough, say less than 4^^^ , we get 



sup|/,(x,OI <2^ 



for all t<t. □ 
The last bound we need is to show is the following local estimate for dtf: 
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Lemma 1.7. Suppose that t is a time such that \pi —pnit)\ < oq for all t < t. For any 

\x — pi\ < ao, the following holds: 

swp\dtf{x,t)\<2 t<t. 

X yt 

Proof. We use a very similar argument to the one of the previous lemma. Write the 
solution as (jl.4p . For the linear part we use again classical estimates (c.f. [1], chapter 
V, Theorem 8.1) that imply 

\dj,\ < Md^zH) . 

For the nonlinear part, consider 

Idthl ^ d^j2 < Gia) f |A„(s)| <4G(a)||MU=o(_i,i)Vt, 



where we have used Lemma [TT3l for the estimate of and (|1.9p . 
Choosing t as in the previous lemma, it follows that 

sup \dtf{x, t)\ < 2Mdk^IzHl for all t < t. 
X yt 



□ 



Remark: Note that the time t depends only on some universal constants, and the time 
the free boundary uses to exit an interval. Unfortunately in the proof of Theorem ll.il 
we will need to choose a smaller time step to < t, which depends on the initial data. 

The above lemma gives us an estimate in terms of the L°°-norm of the initial data. 
Is is also possible to estimate the fx in terms of the L^-norm (which is in our case 
constant in time), at the price of a worse denominator. We will not make use of this 
lemma in this section, but later in the proof of global existence: 

Lemma 1.8. Suppose that i is a time such that \pj — p„(t)| < a/4 for all t < i and that 
sup\ifn)xix,t)\<2Mlhlt±ll^ 0<t<t, \x-pi\<a/i, (1.9) 

X [vty 

then 

<t <i, \x- pi\ < a/4, 

<t <t, \x-pi\ < a/4. 

Proof. The above estimates can be proven following the same steps as in Lemma 11.61 
and II. 71 Instead of using the regularity estimates for the caloric function with the L°°- 
norm of the function, we make use now of the L^-norm. From (Theorem 9, section 



SUp\ifn+l)xix,t)\ < 2' 

X 

SUp\{fn+l)xxix,t)\ < 2 



(Vty 
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2.3.C.) we know that 



sup 

X 

sup 



dh 


< 


ll/IU. 


dx 




d^h 


< 


II/IIli 


dx^ 


t5/2 



where Ii is defined as in Lemma 1 1.61 The rest of the proof follows similarly. 



□ 

In the next corollary we show that given Pn{t) S {pi — ao/2,p/ + ao/2), also the free 
boundary p{t) of the next iteration step is well defined and moreover it is C^^^-Holdev 
continuous. In particular, also p{t) G {pj — ao/2,p/ + ao/2) for t small enough. 

Corollary 1.9. There exists a time to > such that the solution f{-,t) constructed in 
(|1.4[) has a unique zero in [pj ~ aQ,pj + oq), for each fixed time t < to, which we denote 
byp{t). Moreover, p{t) satisfies 

\p{t) -P/| < ai 

where Oi only depends on the initial data as shown in the proof. 

Proof. First, note that CoroUarv 11.51 implies monotonicity of /(•, t) in (p/ — ao,pi + ag) 
which assures the existence of at most one zero, p{t). Next, since f{p{t),t) = for all 
t, then 

On the other hand, CoroUarv l 1 .91 gives a bound for the velocity \ft{x, t)\ < 2 



Wfih 



in the time interval (0,t), and Corollary 11.51 bounds the slope \fx\ from below in an 
interval \x — pi\ < ao in (0,to). Existence results for the ODE in Eq. (|1.10p implies the 
existence of this p{t). 

As a consequence, the interval \p{t) — pi\ is boimded by 



IMI 



L~(-l,l) 



Vi, for < t < to. 



The thesis follows by choosing the constant ai = gMlM^L-Lll^ and the time to as 



rr ao 
/to y. 



(1.11) 
□ 



Now we are ready to show that the map $ defined in Eq. (|1.3p is a contraction in 
the space L°° ((0,to);X), and consequently to define to. At this point we already have 
a specific condition on the time to: it must satisfy all the conditions from the previous 
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results, of which the most restrictive is precisely Eg. p. lip . 







Proof of Theorem ll.lt 

First note that if /„ e L°° {[0,tQ], X) with fn{pn) = satisfies 

PnW e (p/-ao/2,P/ + ao/2), sup i)| < 2^^^^^^^, f„^{p„,t) < 

X V t * 

(1-12) 

then the previous results assure that the next step in the iteration f{x,t) has a well 
defined zero p{t) and satisfies also the same estimates. 

Let /„, g„ e L°°{[0,to],X) be such that fnjpn) = 0, g„(g„) = 0, A„ = -~fn{pn), 
£,n = ~9n{qn), and that both satisfy Eq. ()1.12p . Then we estimate the difference of the 
images by our mapping <!>: for any x G {pi — ao,p/ + ao), 

f{x) - g{x)= I (r{x,Pn - a;t - s) -r{x,pn + a;t - s)) A„ (s) ds 
t 

(r(x-, qn - a;t~ s) ~ T{x, g„ + a;t - s))Cn(s) ds 



t 

[r(a;,p„ - a;t - s) - T{x,Pn +a;t~ s)] (A„(s) - (,n{s)) ds 



+ / £.n{s)[T{x,Pn - a;t - s) -T{x,q„ - a;t - s) 
Jo 

+ T{x,pn + a-t — s) - T{x, qn + a;t — s)] ds 
= ■■ A1+A2. 

It holds 

\Ai\ <snp{t)\Xn{t) - £,n{t)\ / \T{x,Pn - a-t - s) -T{x,pn + a;t - s)\ ds. 

Jo 

We now need to estimate the difference |A„(t) — £,n{t)\ in terms of the quantity ||/n — 
5« li L - ipi-ao ,pj +ao) • Therefore 

\Kit) ~ Ut)\ <\fnMt),t) ~ {t),t)\ 

<\fnx{Pn{t),t) - 9nx{Pn{t),t)\ + \9nx{Pn{t) , t) - 9nxiS.n{t) , t)\ 
<\\fn - 9n\\c^pi-ao.^pi+aa) + \\9nxx\\ {pi -ao,Pi +ao)\Pnit) - gn(0|- 

Next we estimate the term WgnxxllL^^ipi-acpi+ao) with Lemma [1.71 and the C^- 
norm of /„ — with classical estimates for the caloric function. Recall that the esti- 
mates hold in a neighborhood of the free boundary far away from the source and sink. 
Therefore 

|A„(t) - Ut)\ < 'l^"-g"'l^^^^-°°-^^+°°) + 2Mi^zM)|p„(i) _ 5„(i)|, 
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for each fixed time t. On the other hand, by an elementary geometrical argument, it 
holds 

4 

\Pn{t) -qn{t)\ < SUp|/„(x,t) ~ gn{x,t)\ ■ — , 
X A/ 

since the slope is bounded by ^ < |/„^(.T,t)| < 2 



Hence 



|A„(t)-C„(t)| < 



1 + ^ll//IU==(-i,i) 



ll/n 9n\\L°°{pi-ao,Pi+ao) 



(1.13) 



Consequently, the term Ai can be estimated using (|1.13|1 and Lemma 11.31 as 



1^1 1 < 2G(a) 



1 + 



ll/ri 9n\\L^(pi-ao,pi+ao)Vi. 



Similar arguments can be applied in the estimates for A2: the hypothesis (jl.l2p 
gives that 



\A2\<2 
+ 2 



ll/j|iL°°(-l,l) 

Vt 

li//||L°°(-l,l) 

Vt 



\T{x,p„ - a;t - s) - T{x, Qn - a;t ~ s))\ ds 



\r{x,pn + a;t ~ s) — r{x, qn + a;t — s))\ ds 



Next, note that (p„ — a, g„ — a) n {pi — ao,pi + qq) ~ 0. Here we are slightly abusing 
the notation by assuming that Pn < Qn- It holds, for each t < to, that 



\T{x,pn ~ a;t~ s) - r{x, qn- a;t - s))\ ds 



<\p7i-qn\ sup / \Ty{x,y;t " s)\ds 

y&iPn—a.qn — a) JO 

4 

In conclusion we have 

|/(a;) - g{x)\ < cVt\\fn - gn\\L-^(pi-ao,pi+ao), 

where the constant c is given by 



c := G(a) 

Choosing to small enough such that 



/||L°=(-1,1) 



G(a) 



g 



to < 1, 



the proposition is proved. 



□ 
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Remark: We remark that the proof in for short time existence - almost exactly 
the approach of [13] - made use of semi-group methods. The derivation here, while 
admittedly less sophisticated, is certainly more robust and is deliberately tailored to 
the upcoming developments. In particular, we now have tangible criteria under which 
short time existence is purported to break down. Indeed, a parallel derivation for the 
regularized problems allows, at least for the compact case, an immediate proof of global 
existence: With the asset of regularization, flux bounds easily imply various derivative 
bounds which are the central objective of the next sections. 

2 An bound on the flux (and the rescue plan at 
the Neumann boundary) 

In this section we show that, as long as the solution of (P) exists, the flux of the solution 
at the zero set stays bounded. This result will be then used in the next section to provide 
further estimates on the derivatives of the solutions. As a corollary, we will also show 
that the zero set of the solution cannot approach the Neumann boundaries too closely. 

To facilitate matters, we shall, in essence decouple the positive and negative pieces 
of / and, in addition, describe problems of this sort on a larger space which restores 
much of the linearity usually associated with diffusion problems. 

Thus, first, we shall define (P') to be a one-sided version of the system (P) that is 
to say (i) The positive part of / is set, identically, to zero (and no source), (ii) The zero 
p{t) is predetermined. Thus, in principle, in order to recover the system (P), two such 
(P') models can be glued together subject to "additional constraints" on their mutual 
p{t). For the single sided problem, usually to the right of p{t), we will use the notation 
Pp{x,t) for the density, Mp for the total mass, etc. When two such models are to be 
used in tandem, the one on the left (now representing the positive part of /) will be 
denoted by pt- 

Secondly, we may define these sorts of systems - (P') will be sufficient - on a foliated 
space. Let N denote the natural numbers, including zero, and consider [—1,+!]^. We 
shall refer to the individual elements as levels and, denote these, along with various 
associated quantities with a superscript: [—1, +1]*^°^ [—1, +1]'^^\ . . . Let p{t) : [0, T] 
(— 1,+1) denote a continuous function. On each level, we have a copy of p{t) (al- 
ways "located" in the corresponding position) and we consider a sequence of densities 
p^p\x, t), p\^\x, t),..., with the n^^ density supported in \p{t), 1](") C [-1, 1]("'. Ini- 
tially, pp"'(x,0) = for n > while p^\x,0) = pp{x,0). Each of the pp"'(x,i)'s obey 
the diffusion equation with a source to be described below, Dirichlet boundary con- 
ditions at their respective p{t) and Neumann condition at the corresponding x — I. 
Finally, each p^\x,t), n > 1 has a source which is located at its respective p(t) + a 
(where, we remind the reader, a = min{a, ^[1 — p(0]}) ^^'^ strength provided by 
A(")(i) = Vp^"-^'(p(t),i). 

In the context of the foliated model, we may calculate various quantities for the 
original model. Of particular relevance, it is seen that: 

oo 

5]p(")(a;,0 = Pp(a:,t) (2.14) 

n=0 
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and 

/ \{t')dt' :^^n p'^;'\x,t)dx. (2.15) 

To vindicate the above claims, it is first necessary to demonstrate that for t G [0, T'] 
with T' < T, the tail of [—1,+!]^ is sparsely populated. This is elementary and 
follows mutatis mutandis from arguments in the proof of Proposition [2Tl] where detailed 
estimates along these lines are provided. Then, defining 

N 

(2.16) 

n=0 

we have, by monotonicity, 

Pp lim 

N^oc 

satisfying the system (P') which demonstrates Ea. (|2.14p . 
As for Ea. (|2.15p . wc define 

MjrHt)= t pl-Hx,t)dx (2.17) 

Jpit) 

it is seen that J2n Mp^\t) is conserved - and hence identically equal to Mp. Moreover, 
Mp - M^°\t) is given by 

MW(0)-MW(f) = VM^")(t) = f\^^\t')dt' (2.18) 
„>i -^0 

i.e., the mass lost on the first level is exactly that which fluxed up to the higher levels. 
Similarly we have for every level, the identity 

^AfW = f X(^){t')dt' (2.19) 
and summing both sides of Ea. (|2.19p . we obtain Eq. (|2.15[) 

Proposition 2.1. For any T < oo, as long as the solution of (P) exists in [0,T), then 

T 

\{t)dt < oo. 

Proof. As we shall see, the consequences of a finite time divergence violate sensible 
notions of the slow scale for the diffusive transport of substantial material over large 
distances. The pertinent observation is that for St small, there is essentially no diffusion 
over any appreciable distance in the allotted time and in its absence, the "essential" 
supports of the p- and 6-dcnsitics become so widely separated that diffusion cannot 
account for complimentary transports. 

Suppose then that the above display docs not hold. Then for any 5t and K {5t <C 1 
and K ^ 1 to be specified when necessary) there is a ti with ti + 5t < T such that 

pil+St 

/ X{t)dt > K. 
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For simplicity let us reset ti to zero and work with t G [0, The midrange objective 
is to show that under the stipulation of a large flux in a short time, p{t) must head 
towards the boundaries. We consider, for the time being, the one-sided perspective. 

In this context it is noted that the foliated problem in the absence of diffusion is 
particularly simple: mass on each level simply gets displaced, as p{t) sweeps through it, 
to the next level at the position corresponding to p{t) + a. with the density, on each level 
supported to the right of the furthest excursion of p. Thus, in particular if a ^ 1, and 
p{t) docs not enter the right rescue zone (namely x > 1 — 2a) there are at most the order 
of 1/a levels that get occupied in the sweep. To treat the general case, let us consider 
for A G (0, a) small but of order unity the system (P'^) which is defined exactly as the 
system (P') but with a replaced by A. (Also, with corresponding rescue plan, but this 
shall not enter into our considerations.) Let us consider, for identical initial conditions 
and identical p(f;) which remains outside the (Pa) rescue zone, a comparison between the 
(Pa) and the (P') systems. We make two claims both of which are straightforward to 
verify on the basis of an underlying particle model based on non-interacting stochastic 
elements. The first claim is that the is that the total flux in the A-system is not smaller 
than that of the usual system. The second claim is that in two such (P^) systems with 
differing zeros: p{t) < q{t) (with both staying out of the rescue zone) that the flux in 
the p-version is not larger than the flux in the g-version. The proof of these claims 
follow from an obvious variant of Lemma 3.2 in [1]. 

Now let us suppose that p{t) - in a (P) or (P') model - is such that on [0, St], (with 
5t small to be specified later) 

maxp(i) < 1 - 4A. (2.20) 

Then we claim that the flux is limited (e.g., by the order of A"-'^). We use a comparison 
a Pa system with its own q(t). Here g(0) ~ p(0) and then q (almost) immediately jets 
out to a; = 1 — 2 A, then backs to 1 — 4A and stays there till time 5t. Let iVA < 2/ A 
denote the number of boxes in the initial surge which get occupied. Then, by Ea. (|2.15p . 
the initial portion of the flux is not more than N\Mp 

Notice that (still in the g-system with parameter A) all the mass on various levels 
is trapped in the respective regions {x > \ — 2A} with the zero a distance 2A away. 
Thus, in the remaining time, the remaining flux is determined by how much can flux 
across various neighborhoods of size 2 A in time 5t. To this end let us now fix 5t "small 
enough" . To be specific, ((5t)^/^ = 7 A for 7 sufficiently small ensures that the amplitude 
for diffusion across 2A - a unit mass at 1 — 2A and a Dirichlet zero at 1 — 4A is less 
than 

e :^ e(<5t) = e-l/[^*l''^ 

We may also interpret e probabilistically: Starting at the fc'^ level, a fraction (less than) 
£ of the initial mass makes it to the k + T^' level in time St and a fraction (more than) 
1 — e stays behind. Of the former, a fraction (less than) e gets promoted to the k + 2"'^ 
level etc. In short, the upward distribution of diffused mass after time St is dominated 
by the a geometric random variable, X with parameter e and so the total flux coming 
from a unit mass which at i = 0"*" was on the /c*^ level in [1 — 2A, 1] is bounded by 

E(X) = e(l -£) + 2£2(i _ e) + . . . = (2.2I) 
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This latter quantity - independent of k - must be multiplied by the mass that was on 
the k^^ level at t = 0+ and summed. Since the total mass on all levels adds up to Mp 
this provides a bound on the diffusive contribution JXdt that is given by 

i^. (2.22) 

Thus, back in the real problem for fixed A > sufficiently small (i.e., compared with 
the minimum of a and (1 — a) as will be clear below) we may choose 6t small enough 
and take the content of Eq. (|2.22p and add this to our MpN^, and then stipulate the 
sum of these two to be less than one half of the total K. Under these circumstances, 
there must be a first time [St]p < St such that p([(5t]p) > 1 — 2A and moreover, at this 
time [St]p, the total flux is less than ^K. 

With the above in hand, we incorporate the constraints of the full problem - namely 
a pp and a pb both in play with fluxes at p{t) that are supposed to match. The flrst 
implication is that there is another such time [5t]b at which all of pb is in [— 1, — 1 + 2A] 
and moreover, even at the later of these two times, there is more than half the fluxing 
left to be done. Let lb denote the set 

Ib = {te [[St]b, 5t] I p{t) <-l + 4A} (2.23) 

and similarly for Ip. (Here, finally, we choose A small enough so that the intervals of 
size 4A about x ~ ±1 do not intersect) It is clear, from the 6-perspectivc that Jj^ Xdt 

is small. Specifically, an estimate of the form in Eq. (|2.22[ ) applies. Indeed, at the 
earlier time all of pb was confined to the region a; < — 1 + 2 A and has less time, for 
t G than St to achieve diffusion across the gap of length scale 2A. Similarly, from 
the p-perspective the flux during 1^ is small. Since C Ip or vice versa, it is seen 
that there is no time after max{[(5i]p, [St]b} in which the requisite remaining flux can be 
achieved. □ 

As an immediate consequence, we may conclude that - at least in any finite time 
interval - p{t) stays away from a neighborhood of the Neumann boundary. 

Corollary 2.2. Let T < oo denote any time up to which the solution to (P) exists. 
Then there is an e > (depending onT) such that 

p{t) e [-1 + e, 1 - e] forQ<t< T. 

Proof. Suppose not. We may then assume that p(t) approaches —1 as t s- T. In 
particular, there exists a sequence tn ^ T where 

(i) P(tn) = min[t„_t^+5_^] p(t) for some 5„ > and 

(ii) e„ = p{tn) + 1 satisfies e„+i < |;e„. 

This means that the new excursion of p{t) went past the previous excursion of the 
source. 

Let gn-i '■= ^ Xi{t)dt, where Xi{t) is the flux from the first term in the foliation 
model introduced above, with initial data f{x, tn-i). In the context of the frozen model 
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(ignoring diffusion effects), all the mass between p(t„) and p{tn-i) is fluxed through the 
zero set at least once. On the other hand the removed mass via the flux at the zero set 
is then deposited onto the left side by the source term. Therefore we have 

Qn = / pb{x,tn-l)dx > gn~l 

J[p(U).,p(U^^)] 

which would be an obvious violation of Lemma 12.11 

We shall show, more or less, that the diffusion cannot alter this situation in a 
significant way. Indeed we claim that either g„ > gn-i or gn + gn-i > cq, where cq 
is independent of n. Thus, at least one case occurs for infinitely many n, which again 
violates Lemma [2Tl 

Let us consider the situation at time tn~i' a fraction of the mass pb has been fluxed 
and now resides in the interval / = [— 1 + |e,i-i, — 1 + Cn-i]. Let r be the time between 
tn-i and tn where p{t) + 1 = jtn-i- From now on, for simplicity, we denote e := e„_i. 

Note that 

gn > \{t)dt = /(x, U,)dx + /,;_^ -U-l + e/4, t)dt 

> Qn^i + /[_i+£/4 tn-2)dx + -/x(-l + e/4),t)dt 

g-n-l + /[_i+,/4 f{x, tn-.2)dx + A-B, 

where A is the flux of / diffused into / from the region —1 < x < — l + e/4, and B is 
the flux which diffused out from / into x < —1 + e/4. 

Case 1. First observe that if tn — tn-2 < Cie^ then there is not much diffusion 
between [—1 + e/4, —1 + e] over the time interval such that 

B < I f{x, t„-2)dx and g^ > g„_i. 

J[_l+e/4,-l+e] 



Case 2. Next suppose that i„ — t„-2 > Cie^. Our goal in this case is to show that 
either A > B ot gn + gn-i > cq. 

Let a„ := J^l_^_^^2 Pbi^i tn~2)dx . Observe that, by conservation of mass, Mb — an = 
J_^'^'^^^ Pb{x,tn-2)dx. It is then not hard to see that 

A-B>J^ -h^,{-l + e/4,t)dx 

where hix^tn) = an5_i^^/2 + {Mb — a„)(5_i and h solves heat equation in [—1, 1] x 
[in-2,T] with Neumann boundary conditions at \x\ = 1. 
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The effects from .t = 1 are not significant and may be neglected for the time be- 
ing. The condition at x = —1 is, of course, important: the equivalent situation (via 
reflection) is solving heat equation in R x \tn-2, t] with initial data 

h{x,tn) = a„(5„i+e/2 + a„(5_i_e/2 + 2(A4 - a„)(5_i. 

If we assume that a„ < Mb/2, then the flux from the sources at x = — 1 and x = — l + e/2 
regulate each other, yielding 

+ e/4, t) <OioTt> t„. 

Therefore A > B. On the other hand if a„ > Mb/2, then the flux to the a; > — 1 + e/2 
is already significant so that we have + .9n— i 

bigger than cq. The desired result is 

established. 

□ 



3 Global-time existence of the solution 

In the previous section we showed that the problem has unique solution for a small time 
interval tQ. In this section we will show that we can iterate this process to produce the 
unique solution of our problem for global times. 

Let us restart the process at the time to- This will lead to the existence (and 
uniqueness) of a solution in the time interval (to,ii); wc shall continue to iterate the 
process as long as we can. 

The length of the time interval t„+i — t„ in which the contraction can take place is 
a function that depends on 



Suppose there exists a blow-up time < t* < oo: it means that we can no longer 
find a small time interval during which Theorem 11.11 holds . This happens exactly when, 
for any sequence t„ — > t*, one of the following holds: 

(i) limsup„_„„ ||/(-,t„)||^oe(_-^_i) oo or 

(ii) liminf„^oo A(t„) 0, or 

(iii) limsup„_o„ \fxx\{p{tn),tn) oo. 

At the time t* two possible configurations may happen: the limit of the free boundary 
p{t) is an unique point as t —i- 1* , or a puddle of zero forms at t = t* . In both cases we 
will show that the arguments lead to a contradiction, showing the non-existence of such 
at*. In the non-puddle case (section 3.1) a contradiction will be yielded by proving that 
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all relevant norms given in (i)-(iii) are bounded up to < = i*. In the puddle case (section 
3.2) the strategy is a bit different. We first show that (essentially) all derivatives of / 
are bounded up to i = t"^: we then will show that such regularity result is too strong 
to hold at a blow-up time, therefore concluding that the blow-up time t = t* does not 
exist. 

We start showing what happens when the free boundary has an unicj^uc limit. 



3.1 Non-puddle case: limj^t* = 

In this case we can show that, due to classical estimates on caloric functions, none of 
the three breaking factors listed above takes place, therefore yielding a contradiction. 

Let t* the breaking time. When the limit is unique, the source p{t) — a and the sink 
pit) + a also have an unique limit as t — > i*. Moreover, the locations of the source and 
sink at the time t ^ t* are a-away from the point x ~ p{t*). 

This allows to find a spacial neighborhood of p{t*) and a time interval (before t*) in 
which the source and sink never entered in that neighborhood during that time interval. 

Let 

8 

and define Is as a small neighborhood with radius 5 of the point p{t*), i.e., Ig := 
{p{t*) — S,p{t*) + S). Then there exists < ts < t* such that Is x [ts, t*] is far from any 
sink or source. 

In terms of estimates, the fact that sink and source are far away from Is implies 
that we can use classical estimates for the homogeneous heat equation, as in the proofs 
Lemma 11.71 and 11.61 



Let us start the problem (P) with initial datum // := fix^tg). Following the same 
contraction argument for the local existence, we can show that there exists an unique 
solution f{x,t) in the time interval {ts,ts + Ati). Iterating the process another time, 
we can extend the solution to a bigger time interval (t* — J, — 5 + Ati + A<2), with 
At2 > 0. The extension process could continue by iteration. The following is due to 
local estimates for solutions of the heat equation: 

Lemma 3.1. Let f{x, t) solve (P) in the time interval (ts, t* — S + A). Then 
\f(^)ix,t* - S + A)\ < C^^^^^''''^' znls. 



where f'--'^^ denotes the k-th spatial derivative of f. 

Due to above lemma, the mass conservation, and the iteration process, we obtain 
+ < (A/p + A/fc)(r)-('=+i/2) for < r < t*-ts. 
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In particular X{t) is bounded, and therefore the / is uniformly bounded up to i = t* . 
Therefore case (i) and (iii) are eliminated. 



Moreover ft is also uniformly bounded near t = t* in Ig, and thus f{x,t) uni- 
formly converges to a C°° function F{x) as t ^ t*. Let us call this function f{x,t*). 
Then f{x,t) solves the heat equation in Ig x [t^jt*]. Below we show that F{x) is non- 
degenerate at X = p(t), thus eliminating (ii). First let us show the analyticity of F in 

Lemma 3.2. For every tg < t < t* fixed, the Junction f{-,t) is analytic in Ig. Moreover 
the function f(-,t*) is positive to the left of p{t*), and negative to the right of p{t*). 

Proof. We write 

r + l 

f{x,t)^ J ^ r{x,x';t)fi{x')dx' 

+ f [r{x,ip{t')-a);t-t')-r{x,{p{t') + a);t-t')]X{t')dt' 
Jo 

The part of the solution coming from the linear term is analytic (see [5] Section 
2.3.c). For the nonlinear part, it is an easy computation to show that the function / 
extended to the complex plane f{x) — > f{z) with z = x + iy is complex diffcrentiablc 
i.e., satisfies the equation 

dx dy 

As for t = t* , note that the formula still holds for F{x) ~ f{x, t*). 

For the second claim observe that f{x,t) solves the heat equation in the domain 
Is X [tijt*] with continuous zero set p(t), and with initial data f{x,ts) which was 
positive to the left ofp{ts) and negative to the right otp{ts). Furthermore on the lateral 
boundary of Ig, f{x, t) > on the right and f{x, t) < on the left. Therefore the claim 
follows directly from the strong maximum principle for solutions of heat equation. □ 



We try now to take care of the possibility [ii] . Assume that 

f,{p{t*),t*) = 0. 

The following holds 

Lemma 3.3. Let t* the first time at which the gradient at the free boundary degenerates, 
i.e fx{p{t*),t*) = 0. Then, also 

/,,(p(r),r) = o. 



Proof. The proof follows from a very simple observation: Suppose that fxx{p{t*),t*) > 
0. This fact cannot coexists with the properties of the hmction /: / is positive for x < 
p(t*) and negative for x > p{*) with purportedly vanishing derivative and very smooth 
across the free boundary (in fact / is spatially analytical in space in a neighborhood of 
the free boundary, as we have shown). □ 
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The question is now how to define the velocity of the free boundary at the time t* . 
Note in fact that as long as fx{p{t),t) stays away from zero, the speed oi p{t) is given 

by 

which, as we have just shown, is an indefinite form if the first derivative degenerates at 
x=p{t). 

Then, in order to find define the speed of p{t) at this critical time t*, let us take the 
Taylor expansion of the equation f{p + Ap, t + At) = 0. It holds 

0=f{p + Ap,t + At) = f{p,t) + fx{p,t)Ap + ft{p,t)At 
+ ^[f.^Ap^ + 2f,tApAt + fttAt^] 

+ ^[/xx.Ap3 + 2/,,tAp2At + 2f,uApAt^ + ftuAt^] + - 

At the time t*, from the Taylor expansion we get that 

f 

XXX 

We would like now do determine if this speed in finite or infinite: from the Hopf Lemma 
we know that in a parabolic problem, as long as the free boundary has finite speed, the 
derivative at the boundary is strictly positive (or negative) for all time. This implies 
that in our case, at the time t* the speed of p is not finite. Consequently the value 
fxxx{p{t*),t*) has to be zero. Now we can repeat Lemma 13.31 for the derivatives fxxx 
and fxxxx and prove that also 

fxxxx{p{t*),t*)=0. 

Again using the Taylor expansion we would then define 

f 

p = ~^f^, 

•J J xxxxx 

etc. Repeating Lemma 13.31 iterativelv. in is seen that all the derivatives of the function 
/ at the time t* vanish. 

Since all the derivatives of f{x,t*) are equal to zero at p{t*) from the previous 
lemma, analyticity implies that f{x^t*) is zero in the entire neighborhood Is as well as 
the point p(t*). 

This is a contradiction since we have showed in Lemma 13.21 that any solution at the 
time t* is strictly positive for x < p{t*) and negative for x > p{t*). 

Therefore the gradient of the function at the free boundary p{t) stays strictly nega- 
tive for any time t > 0. 

This eliminates case [ii] and we are done if there are no puddles. 



3.2 Puddle case: C := liminft^t* <1Z := limsup^^^* 

Here the situation requires more careful analysis: first we claim that if a puddle forms 
at i = then all spatial derivatives of / are bounded in the whole domain, at least 
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along a sequence of times converging to t*. (For precise statement see Proposition 
13.41 below). This result, along with an analyticity argument, will immediately yield a 
contradiction. It is remarked that in the regularized case, such results follow, essentially 
from quadrature. However, these estimates deteriorate as the rcgularization is removed 
so wc shall not pursue this venue. 

In the case of a puddle, it is clear that 

J ^ as t T (^*)- 

(This is because the zero travels between the endpoints of the puddle so fast that there 
is not enough time for mass to diffuse into or flux through the zero set as t t* .) 

We begin by choosing to close enough to t* such that 
(BH) \f\{x, t)dt < ^ min(M6, Mp). 



3.2.1 Derivative bounds 
Proposition 3.4. Suppose 

liminfp(i) ~ C <TZ ~ limsupp(t). 

Then there exists a constant C (m) such if P > is sufficiently small then the following 
holds: for any sequence tn — > t* such that p{tn) — > Ti- 

limsup|/('")(x,i„)| < C(m) for any x [U - [3,11 + [3]. 

n — >oo 

A parallel result holds for C 



Proof. Since p(tn) Ti, the number j3 is chosen so that a{TZ) > 3/3. This means that 
for all n sufficiently large, 

p{t„) + a{p{tn)) >n + 2l3 (3.1) 

i.e., during the times of specific interest the sink is well outside the interval under 
consideration. Moreover at some point t in [to,t*], p{t) — a{p{t)) < TZ — 2/3 for t > t. 
Thus, for all times after i, for any xq G [R — (3,TZ + (3], the source only represents a 
distant agitation. With this in mind, let us reset to such that all of the above holds for 
all n and/or for all times greater than Iq. Moreover, due to an up and coming plethora 
of indices, we might as well define i„ =: i*, with (only) the property of Ea. (|3.ip to be 
reserved for later. 

It will prove convenient to work with with the auxiliary variable which measures the 
time remaining: 

9 -.^ 0{t) = t* - t (3.2) 
for t € [to,i'^]. Let B{6) denote a slowly diverging function e.g., to be explicit 

B=^\\og9\^/^ (3.3) 
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will be adequate with Bq chosen for the purposes at hand. For G [TZ — P,Tl + P] let 
us estimate magnitude of the to"^ derivative of /(xo, i*). 

We denote by Xq := XQ — a which locates p{t) when the source contributes to /(xq, •) 
and its derivatives. If we denote by C{t) the distance between p{t) to Xq then by the 
Green's function formula 

\f-'\xo,t*)\<J^ r^e-^(*)^/«dt (3.4) 

for some L = L{m). We start by defining the set He [to, i'^] via 

M = {t\C{t) < B{e)Ve}. (3.5) 

Our first claim is that H is the only important set for the ostensible development of 
singularities in /(xq, •). Indeed 

f ^e-^WV»dt< / ^e-^'dt< f* \{t)e^«-^dt (3.6) 

which converges and is small independent of t'^ < t* for Bq large enough. Henceforth 
we may focus on events that take place when t gM where, for all intents and purposes, 
there is no help from the exponential factors. (But, on the positive side, H is disjoint 
from the tail end of [io,t*]) 

To aid with our objectives, it will be convenient to divide [to, t"^] into disjoint regions 
that are of equal size on a logarithmic scale: Let 3> 1 denote a sufficiently large 
number the precise (minimum) value of which will be determined in what is to follow. 
Roughly speaking, we wish the k^^ region to be of size H^^ of the fc — 1^'. Specifically, 
we may proceed as follows: The fc"^ region will be denoted hy gk, fc = 0, 1, . . . and the 
size, \gk\, will satisfy = i?-'=|.go|. (Thus |go| = ^{t* ~ to)). Therefore 

go = [to,to + |go|), ••• ,5fe = [to + l5o| + --- + |,9fc-i|,to + |go| + --- + l5fe|]- (3-7) 

Note that the value of 9 at the right end of gk is a constant (that is very near one) times 
the size of gk+i- 

We shall also need two spatial regions: A2, A3 which are given by A2 = {X^ + b, Xo + 
2b) A3 = {Xq + 2b, Xq + 36) with b chosen small enough so that these sets are inside 
[£, TZ] lying well (on the scale of b) to the left of TZ. Moreover b is large enough (or to 
late enough) so that b is much larger than B{9{to))y^6{to). We shall define epochs that 
are punctuated by certain exits from and entrances to the region A3. The beginning 
of an epoch, denoted by T,„i,i is when p{t) enters A2 from A3 and will not revisit A3 
before first having touched the appropriate B\^ neighborhood of Xq- Let r^un be the 
first time during this epoch where p(t) visits Xq + B^/O. The time Tmax is when p{t) 
leaves the ostensibly smaller B^/O neighborhood and does not touch these neighborhood 
types till another visit to A3. To be specific, Tmax will be the moment of this departure, 
so p(Tmax) ~ Xq + [Sv^] [rmax ■ There is a third time, namely r,„in when p{t) actually 
enters the appropriate B^/O neighborhood but this time does not play a particularly 
major role. 

Our first claim is that if epochs are localized e.g., to a single gk then their contri- 
bution to [[/'^"'^(xo, •)[[ is tractable. Indeed, assume for simplicity that times t[ji|,j < 
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'^min < ■''mix < ' ' ' < T^^^^ < '^min < "^mix are the Only punctuation marks in gk- Then 

where 9k+i := — (to + •■• + \gk\) ~ l5fe+i| is the time remaining by the end of gt- 

Next we show that the integral to be done is actually of the order of exp[—l/9k]. 
This will be facilitated by the perspective of pp: Note that at t = T^niin^ ^^"^ density 
Pp{x,T^^-^^) vanishes for x < Xq + 2b and has some positive profile for x > Xq + 2b. 
Starting from this profile, we are supposed to compute the flux of pp through p{t) while 
t G M. It is not hard to see that this is less than the total flux through (the vicinity 
of) Xq with Dirichlet boundary conditions at x = Xq + B^fO, and an initial profile of a 
delta mass at Xq + 2b. The strength of the 5-mass might initially taken to be the total 
of Mp. But to account for the possibility of refiux from the p-source (sink) - which is 
much further away, let us add to this the quantity 

A := / \dt < oo 





which is all that ever has fluxed and ever will flux. 
This leads to the estimate 

Xit)dt < C2C-'=^/^'= (3.9) 

where C2 < oo and C2 > arc constants which do not depend on fc, or the end time . 

The key issue is therefore to show that the epochs do not extend over many g- 
scales. Indeed, it is remarked, an extension of the above reasoning tracking a single 
epoch through scales k < k < K would yield a bound of the form 

gle--/«» (3.10) 

replacing the corresponding right hand side of Eg. p.Sp . And, as unlikely as it may 
seem, if i^T ^ k this could be large. Wc further remark that some positive powers of 
C{t) = \Xq — p{t)\ originating from various places are available (inside the integrand) 
for the estimate. But even in the best case scenario - namely the estimate for the norm 
of /(x'o, t) itself ~ there is still a logarithmic divergence (oc K — k) weighing in against 
the exponential prefactor. Finally, it might be technically noted that since t* < t* it 
must be the case that all but a finite number of the gk are devoid of epochs. However 
the technicality is of no practical significance since we seek bounds that are uniform in 
< t*. We turn to the task at hand. 

Consider an epoch defined by the times fmin and Tmax- Let A3. 2 denote the left and 
middle third of the A3 region namely 

A3,2^{x\Xo + 2b<x<Xo + h} (3.11) 

and let us estimate the accumulation of mass in the region A3 2 at time t = t^^x- Note 
that in addition to the flux from the outside of the region, there is the initial mass that 
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was in A3 at time t = fvnin- This is actually of negative utility. We cannot rely on 
it staying in the region and we do not wish to account for where it might go during 
Tmin < t < Tmax- Howcver, in accord with condition BH, this could only account for 
10% of Mp. Then by mass conservation there is a significant portion of Mp outside of 
to the right of TZ. 

Thus, we shall rely on this material that diffusing in from the right. Here, it is 
worthwhile to recollect that the underlying condition of the epoch is that p{t) stay out 
of A3. So, placing the guaranteed fraction of Mp at the extreme right - x = +1 - and 
placing Dirichlct boundary conditions at x ~ Xq + 2b wc obtain the estimate 

Q3.2= f Pp(a;,w)rfa;>C3,2C-^^'^/[^(^— (3.12) 
JA3.2 

where C3_2 > and C3 2 < 00 are constants independent of the parameters of the epoch. 

Now just about all of Q3^2 and more will be swept into p{t) between time Tmax and 
the next time that p{t) crosses all the way to the right side of the region A3. This may 
happen immediately but it might incur the passage of other possible (long) epochs. 
However, from the stipulations about the time sequence (t„) it is inevitable that this 
will happen before t = . 

Our final substantive claim is that Q3,2 is essentially dominated by the flux through 
the left boundary of A2 in the time interval [rmaxji'^] by a delta-source of strength 
Mp + Mt + A placed at Xq + B^/9\r„^^^ at time Tmax- The implication provides an upper 
bound which we state as a separate lemma: 

Lemma 3.5. Let Q3,2 be as described. Then 

Q3.2 < C4e-^*/«(^''-) (3.13) 

where C4 and C4 do not depend on the parameters of the epoch nor on the time t"^ . 

On the basis of Lemma 13.51 and Eq. (|3.12[) , wc may conclude that 

) > Ci0{ ) (3.14) 

with ci > independent of the parameters of the epoch and the time t"^ . Thus consider 
H to be sufficiently large {H > 1 /ci) and let us double cover [to, t"^] with two overlapping 
partitions of the general type described - both with parameter H - in such a way that 
(on the logarithmic scale) the midpoints of one partition form the endpoints of the 
other. 

Then according to Eq. (|3.14[) . all epochs are caught in individual elements of one 
or the other partition. This enables us to sum (twice) the analogue of Eq. (|3.10p with 
matching fc's all the way to fc = 00 which demonstrates the desired result. □ 

Proof of Lemma 13.51 

Consider the situation at t = Tmax where that there is non-zero pb for x < Xq and 
non-zero pp for x > Xq. The initial densities (and ensuing fluxes) are, evidently such 
that soon enough, p{t) will exit from the neighborhood of Xq and cross all the way into 
A3 . Later - maybe much later on the log-scale - it will cross all the way through A3 

Now consider p{x, t) defined in the domain T,^ := [—1, XQ + 3b] x [rmax, t"^] as follows: 
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(a) {dt - A)p = X{t)6^^p(^t)-a ; 

(b) 9,p(-l,t) = 0; 

(c) p{Xo + 3b,t) = ; 

(d) p(x ) + ~ Pp{x, i)X(Xo+2fc,Xo + §6) i^) ; 

where (f is the mirror image of Pp(-, Tmax) in A3 reflected about x = + b, i.e., 
V3(a;) := Pp(2(Xo + b) - T,nax)X(Xo-6,Xo) (2;)- 

We denote a(<) := {x : a(x,t) = 0} and denote by be the first time when a{t) 
touches Xo + 36. There are two observations: 

• First, by maximum principle for caloric functions it follows that p < p in S"^, and 
thus 

p{t) < a{t) for w <t<ti (3.15) 
This, of course, places < t^. 

• Secondly, observe that p{x,t) > p{2{Xo + b) — x,t) in E, again by the maximum 
principle for caloric functions. In particular, 

e [^0 + ^^0 + 36] =A2U A3; w < i < (3.16) 

Since J p{x,Tniiix)dx = Q2.3 and J p{x,t^)dx = 0, t = t^, it follows that 

Q2,3 = y {-p.)ia{t'),t')dt' + J p,{Xo + 3b,t')dt'. (3.17) 

Moreover since we know at € A3 U A2 for all t S [rmax, by obvious dominance 
the second term in Eq. (|3.17p is bounded by the flux of h{x,t) through x = Xq + 36, 
where h{x,t) solves the heat equation in the region to the right of Xq + b during the 
times [Tmax, with initial data the same as p at t = ^max and Dirichlet conditions both 
sides of A2 U A3. Since we may write that the fraction of (53,2 that is lost via 

the right boundary is less than C'se"'^^/^^'^"'"'^ where the C5 and C5 only depends on b. 
(So that the total which is lost is no more than Q2,3C5e~'^^^^^'^'"'^''\ As far as we are 
concerned it is sufficient that this is less than half of Q2,3-) 

It remains to estimate the first term in Ea. (j3.17p . For this purpose, let us consider 
g{x,t) solving the same problem (a) — (d) in S as p, except that g{x,t) has Dirichlet 
boundary conditions at a; = Xq + 6. By caloric inequalities the flux of p through a{t) 
(where a{t) > 6) in the interval [Tmaxi^"''] is, from the perspective of the left, less than 
the flux of g through the line x = Xq + b. Finally we may further modify g to a g 
which has the Dirichlet condition at a; = Xq + b has all conceivably available mass ~ 
namely Mp + Mt - is placed at Xq. Moreover, the source term is placed as far forward 
as possible and allowed (more than) all of its available flux as soon as possible. This 
amounts to adding A to the Mp & Mb which are in the mass at a; = Xq Thence 



24 



where C4 and C4 only depends on b, Mp, etc. □ 

Now may now we finish the proof of Tlreorem 10.11 
Corollary 3.6. There does not exist a finite blow-up time t* . 

Proof. We already eliminated the non-puddle case in Section 3.1, so let us discuss the 
puddle case. Let £ and TZ be as above. Let t„ — > t* be such that p{tn) TZ- In 
particular, choose n large enough such that \p{tn) — P\ < e ^ |a. Call fn{x) = f{x, tn). 
The function fn(x) has a unique zero at p{tn)- In Lemma 13.21 it was proven that 
(for an sufficiently large) fn is an analytic function in any subset of the open interval 
{p{tn) — ^a, p{tn) + ^a) . This imphes that fn{x) is analytic in the interval {TZ — /3,TZ + f3) 
for large enough n. 

Moreover due to Proposition 13.41 the sequence {/„} is uniformly bounded and uni- 
formly Lipschitz continuous in the 2/3 neighborhood of TZ. Hence Ascoli-Arzela Theorem 
ensures the existence of a subsequence such that /„j. — > 4'{x) as i„ t* uniformly 
in (7^-/3,7^-K/3). 

Now, since the sequence of analytic functions f„^ converges uniformly, the limiting 
function (f> is also analytic in {TZ — P,7Z + P). Now let us observe the profile of near 
TZ. From the choice of the sequence tn, it is clear that 4>{x) is positive to the right of 
X = TZ. (Indeed for any e > 0, for sufficiently large n depending on e the function f{x, t) 
solves the heat equation with source term in S,: := [—TZ + e, +1] x [tn, t*) with positive 
boundary data at TZ + e and t = tn and Neumann boundary data at x = 1 . Therefore 
f{x,t) stays strictly positive in S^, staying uniformly away from zero as t t*. ) 
On the other hand cf) is identically zero to the left of TZ. These two facts interdict the 
possibility of analyticity for the function 0. □ 
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